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Effect  of  Multipath  on  the  Height-Finding 
Capabilites  of  a Fixed-Reflector  Radar  System 

Part  l:Analy$is 


1.  INTRODUCTION 

In  this  report,  we  present  the  analytical  methods  for  determining  the  error 
(due  to  multipath)  in  the  altitude  measuring  capability  of  a reflector  radar  system 
operating  in  a monopulse  mode.  In  particular,  we  will  study  the  case  of  a fixed 
reflector  antenna  on  which  there  are  two  separate  horns  in  the  focal  region,  with 
the  horn  output  processed  (see  Figure  1)  in  such  a fashion  as  to  give  altitude 


REFLECTOR 


Figure  1.  Motiopulse  Operation  of  the  Deflector  Antenna 


(Deceived  for  publication  18  May  1978) 
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information.  The  analysis  we  will  present  here  is  valid  for  an  arbitrary  reflector 
located  at  an  arbitrary  height  above  the  earth  and  for  an  arbitrary  location  of  the 
two  feed  horns  in  the  focal  plane.  In  part  II  of  this  study,  we  will  apply  the  results 
developed  here  to  specific  radar  systems. 


2.  ELEMENT  PATTERN  (NO  EARTH  - THAT  IS,  VACUUM) 

In  the  Physical  Optics  Approximation  (reflector  much  larger  than  wavelength) 
the  magnetic  field  scattered  by  a perfectly  reflecting  surface  S is  (see  Silver1) 


S 


dS(n  X Hj)  XV 


(1) 


where  R is  the  distance  from  a point  in  dS  on  the  surface  S to  the  observation  point, 
n is  the  unit  normal  to  the  reflecting  surface,  and  H is  the  incident  magnetic  field 
from  the  primary  source  (horn  antenna).  Consider  the  system  in  Figure  2.  The 
reflector  surface  is  described  by  the  arbitrary  function 

z s f(x,  y)  . (2) 


The  surface  Integral  in  Eq.  (1)  can  now  be  written  in  terms  of  the  projection  of  the 
reflector  onto  some  plane,  as  shown  in  Figure  2.  If  we  denote  the  projected  aroa 
as  S0,  then  (sinco2  dS  * dx  dy  {1  + (0f/dx>2  + (df/Oy)2) l/2) 


-a  ■ [i  * (-|-)a ♦ ($]  , 


(3) 


1.  Silver,  S.  (1955)  Microwave  Antenna  Theory  and  Design.  Dover  Publications. 
Inc.,  New  York, 


2.  Kaplan,  W,  (1953)  Advanced  Calculus.  Addtson-Wesicy,  Publishing  Co. , Inc., 
Cambridge,  Massachusetts. 


Figure  2.  Assumed  Reflector  Geometry.  Azimuthal  angle  ^ of  observer  is  not 
shown 
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and  x,  y and  z are  unit  vectors  along  the  x,  y,  z axes.  Now 


a --jr*  (x  sin  0 cos  d + y sin  0 sin  d + * cos  0}  . 
o 

• exp  {-  ik{R  - x sin  0 cos.d  - y sin  0 sin  d - z cos  0)) 


(5) 


in  the  Fraunhofer  Zone.  If  we  use  (4)  and  (5)  in  (3)  we  get 


Ss  • w //  d*  d4(J  - - Si] 

O Q ' ' J 


X [x  sin  9 cos  <t>  + y sin  6 sin  <j>  + z cos  9]  (6 

* exp  {-  ikRQ  + ik(x  sin  9 cos  <£  + y sin  6 sin  + z cos  0)}  . 

Note  that  in  (6),  the  quantity  z = f(x,  y).  From  practical  reflectors,  quite  often 

2 4 

z = f(x,  y)  can  be  expressed  as  z = o(x)  + 0(x)y  + y(x)y  . The  problem  now  is  to 
calculate  the  incident  magnetic  field  Hj  due  to  the  horn  antenna.  For  the  horn 
antenna  representation,  consider  Figure  3.  We  characterize  the  horn  radiation 
pattern  by  the  electric  field  distribution  (horizontal  polarization) 


Si » hmo)«K) 


where  i t>Q  is  the  unit  vector  shown  in  Figure  3 and  can  be  written  as 


<&o  = y*  cos  r)0  - z1  10  (3) 

and  y‘  and  z'  are  unit  vectors  along  y'  and  z'\  Also  s0  is  shown  on  Figure  4 and  is 
given  by 

2 2 2 

a0  ■ l(x  - x/  * y*  * (a  - isor|  . <0> 

From  the  geometry  in  Figure  4,  it  can  be  seen  that 


C>0  « tan 


% * tun 


t *0»^»  *0 ) 

b'lgure  4.  Crows  ".Sectional  View  of  Hoflector-Horn  Geometry 

where  tan  _ is  of  modulo  . From  Cq.  (7),  we  see  that  wo  are  assuming 

Kj  is  horizontally  polarized  ami  ran  be  represented  as  the  product  of  the  arbitrary 
functions  (that  is,  they  can  be  anything  we  want  them  to  be)  1>’(0())  and  g(»j  ).  Note 
that  a(),  d,(  and  all  depend  on  x and  y,  as  is  evident  from  Kqs,  (0)  to  (11). 

We  nest  have  to  calculate  the  incident  magnetic  field.  This  can  be  obtained 
from  (7)  via 
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(12) 


^0  =* 


sin  cos  riQ  sin  rjQ 


cos  ^ sin  a + cos  p sin  0 cos^  n 
o ro  o _o 


p = - tan 
Ko 


‘ £) 


If  we  now  use  Eq.  (16)  in  (14)  and  then  use  that  result  in  (6),  we  get 


Ik  Y.  e 
H _ o 

-a  ‘ 2»R 


rr  f(0  ) § (n ) 

JJ  dx  dy  HAox+Boy  + Co£} 


-ik(s  -xsin0cosi£-ysin0sin$-zco30) 
% ^ 


where 


Ao  - COS  0 (%  + Y0-g)  - sin  o sin  * (oQ|y  - 0of£)  . 

Bo  = cos0  fc0  + Y0|r)  +8taOco3*(“o-§’0ow)  • 

CD  - - sin  0 sin  * (po  + *yQ  f^r)  - sin  9 cos  4 (»Q  + 70f~) 


The  radiated  power  pattern  is  proportional  to 


s7t| 


• «o^8  X V *^“3  * V XV 


“ Y^»a  ’ V 


where  n is  a unit  vector  in  the  direction  of  the  observation  point. 


•>:  ’>"v 
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■v’^nseniai^’ot, 


If  we  use  Eq.  (21)  in  (25),  we  get 


47r  Po(0,  4) 


F(ip  )g(n  )A  -ikp 
dx  dy £— — - 0 ° e *o 


where 


/*/*  WJg(n  )B 

JJ  dx  dy 2— — £ S 


ff  axdyf^!i^e'ik«>o  2 i 


p0  s aQ  - x sin  fl  cos  - y sin  0 sin  0 - 3 cos  6 


3.  multipath  effects 

When  the  reflector  in  Figure  1 is  located  above  the  earth  the  field  «t  p <„ 
consist  not  only  of  the  direct  component  in  Eq>  (21),  but  also com 
"»”*  acaiicred  from  ihe  earth.  as  „ indicated  in  Flenre  5 Ftaro  5 IT 
do«od  ,ho  transmitter  altiinde  as  d.  u»  airplane  ai.iinde  as  h,  and  the 

*f  ;:  r~:  ~ 

in  order  to  ealoniate  the  molt, path  field.  „ „ „Moaaar,  M 

r^rrrrcrr  — — - 

(a)  h > d 
‘ i B I + P cos 

v * ' (20) 
3<  Dg^5)  WaveSx  Dover  ^lca(i, , 


rsrl  + r2  • 

For  the  dish  higher  than  the  airplane,  we  have 


(b)  d > h 


D = 


1 + 


4S^S2T 


sd  - spa  + t)J 


-1/2 


(33) 


Sl  = 


(2aed) 


TJl 


S2  = 


(2aeh) 


172 


/hN1/2 

T = (l) 


rl  = r " r2 


Now  from  Figure  5,  we  see  that 
r„ 


vt 


e 


(34a) 


1 " a. 


(34b) 


R2  = [<ae  + d)2  + a2  - 2ae(ae  + d)  cos  63! 


<3  3a) 


»2^-2 


^1  * ^ae  + + ae  2ae^ae  + ^ C03 


1/2 


(35b) 


Also 


0"  o a in 


2R2(ae  + d) 


6 . 


(36) 


We  can  next  write  the  incident  magnetic  field  at  the  reflection  point  Q (see  Figure  3) 
as 


H(Q) 


ikYQe 


“lkR„ 


H, 


// 


F(0  ) g(n  ) -ik(a  +xain0"-zcos0") 
dx  dy  - — 2-e  0 


‘ (A”  i + + €;;£} 


(37) 


>~rxn  W 
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1 

1 

1 

i 

I = Hxi  + Hyy  + Hzz  , (37) 

1 

| where 

| Ai,'“sS"(%ti'o-M-)  • ««> 

i 

| 

By  cos  0"  (Po+,o|I.).stae„(„o|I.PoiI)  , (3I)) 

I 

cy  ■ sta  e..  (oo  ♦ ,o  §J)  . (40) 

Now  we  transform  the  fields  to  the  prime  coordinate  system,  which  is  shown  in 
j Figure  5.  We  have  the  unit  vector  relationships 

j 

x = x1  cos  (02  - 6)  - z1  sin  <02  - 6)  , (41) 

! and 

I 

| 

z = z<  cos  (02  - a)  + x'  sin  (02  - 6)  . (42) 

Also,  the  unit  vectors  and  n^  directed  along  the  incident  and  reflected  wave 
vector  are 

o z*  cos  <J>2  - x*  sin  02  , 

i 

= z1  cos  02  + x'  sin  02  . (43) 

Therefore,  in  the  primed  coordinate  system 

j 

| Ht(P>  » Hx  (x*  cos  (02  - 6)  - z'  sin  (fi2  - 6)}  + I^y' 

+ lf2  {z*  cos  (02  - a)  + S'  sin  (0g  - a)) 

(44) 

* x'  (Hx  cos  (02  - A)  + Ha  sin  (02  - a))  + y'  Hy 
! ♦ S‘SHa  cos  (®2  - 6)  - sin  (02  ~ * 
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and  the  incident  E field  is 


Y0  E^(P)  = x'  H cos  i>2  - y'  [Hx  cos  (^2  + ®2  ' 6)  + Hz  sin  + S2  ‘ 6)J 


(45) 


+ z'  Hy  sin  02  • 


Now  at  the  reflection  point  we  have  that*  (for  perfect  conductivity  or  very  small 
grazing  angles) 

(erV  ■ - ®iV 

<Vx,  * «¥,. 

(Er)x,  means  the  x1  component  of  ER.  Therefore 
(,2' 

- z>  H sin 


Y 

E„  = x<  H cos  0,  +y'{Hv  cos  <(L  + e„  - 6)  + H sin  (0„  + 6 - 6)  } 

Ip<Md“r  y -s  x 2 2 Z 2 2 


(46) 


where  |o(^2)|  ~ 1 is  the  reflection  coefficient  at  Q,  and  D is  the  spherical  disper- 
sion factor  defined  in  Eqs.  (29)  and  (33).  Now,  the  reflected  magnetic  field  is 

&tW£r>  • 


or 


H 


Jn  « 

jg5-  * x'(-  Hx  cos  ^2  cos  (V»2  + ®2  - 6)  - Ha  cos  02  alw  W2  + ^ 


♦y*  H. 


(47) 


+ z'[Hx  aln  <i> j cos  (ti2  + ®2  “ + Hz  aln  ^2  8ln  ^2  + ®2  " 


we  will  assume  that  the  reiiectlnn  is  specular.  This  Ignores  the  effect  of  diffuse 
reflection,  and  is  valid  provided  the  height  a of  the  surface  bumps  are  such  that 

a < 3 where  X » signal  wavelength. 

2 


•4. 


s' 

:r 

1 

? 

| We  now  use  Eq.  (42)  to  express  HR  In  the  (x,  y,  z)  coordinate  system 

i 

i 

\ 

\ H 

1 = - x cos  W2  - «2  + cos  <*2  + fl2  - 6)  + Hz  sin  (*2  + 02  - 6)) 

| 

j 

1 +y  h 

1 y 

« 

(48) 

T 

j 

| + z sin  <*2  - 02  + 6)IH%  cos  (^2  + 02  - 6)  * Hz  sin  <02  + #2  - 6)) 

| 

• 

f 

Ij 

\ 

1 

1 

| Therefore,  the  multipath  magnetic  field  at  the  airplane  location  is 

f 

1 

J 

'i 

! 

{ -IkR, 

I Ik Y De  3 ff  FW- )g(ri  ) -ik(s  +xaine"-zcosd") 

Sr-  °r0  li**  s0  6 

So 

(49) 

i 

i 

i 

J 

] 

\ 

i 

. {A^'x  + B^'y  + C”'^  , 

i 

$ 

i 

j where 

! 

j 

r3  = rx  +r2 

i 

} 

\ 

and 

V 

s 

Ao"  B “ (“o  + volx)  008  ^2  * °2  + 6)  003  ^2  + S2  " 6 “ 0M)  » 

(50) 

t 

J 

■ j ’ coa  fl”  (*o  + To  ly)  ’ 8ln  ®"  (tto  § “ 0o  li)  * 

• \ 

(51) 

i 

.1 

\ 

i 

! Co"  ■ («o  + Yo  lx)  ain  ' °2  + 6)  coa  <^a  + e2  ‘ 6 “ 0n)  • 

(52) 

\ 

i 

Now,  from  Figure  5 it  is  clear  that 

“a *(5  + *2).*  • 

.$ 

or 

j 

^2  ♦ °z  “ 6 " ®"  " 0 • 

(53) 

Therefore 


Am  = - (oo  + yQ  g)  cos  - 02  + 6)  , 

Bo"  - 003  0"  (Po  + >o  §)  - s*  «"  (%  § * 0O  8i 


Co ' = (Qo  + ll)  ain  <*2  ’ ®2  + 6)  • 


We  now  obtain  the  total  magnetic  field  at  P by  adding  (21)  and  (49) 


-TOTAL(P) 


-ikR 

ikY  e 0 


ff  dxdy 


F(0  )g(n  ) -Iks, 


* pto6**1  ^'-.ve‘l54] 

+ y [B0®5  1 + Bo"  ve  ^4] 

♦*  [Co®15  l + Co’ ve’?4]|  * 


where 


Cj  = k(x  sin  6 + z cos  6)  , 

$4  ■ Wjc  sin  0"  - a cos  0"  + R - a ) 

3 0 • 

v-!p<02)|d^d  . 


•wHa-a**-** MfwWfraWtM  I-  -a  .-. . , 


Therefore,  the  power  pattern  is  (including  multipath) 


4*PM0<e> 


rr  f(0  )g(« > -ucs  r -nj 

JJ  dxdy  -S--  — e ° [Aq  e ^A-ve  4J 

So 


FttUgGlJ  -Iks 

dx  dy 2. 1 e 

3o 


0 [b„  «l5‘ + 


B'"  ve 
o 


F(tf-  )g(rjJ  -iks 

dx  dy 2- — e 

3o 


o r ih  2 

[coe  +CJ)"ve  4J  . 


We  now  relate  0 and  h.  From  Figure  5,  we  first  calculate  that 
R2  3 *ae  + d*2  + <ae  + h>2  " 2<ae  + d>(ae  + h)  cos  ©T  , 


rl  + r2  r 

8 1’0 t'-hr1-*: 

e e 


Also,  from  Figui'e  5 we  can  relate  9 and  h via 

<ae  + h)2  - R2  + (ae  + d)2  - 2Ro(a0  + d)  cos  (»/2  -6  + 9)  . 


Therefore 


. (a  + h)2  - (a„  + d)2  - R2 
0(h)  = sin  ■‘2fto(ae  + d)  + 6 * 


Therefore,  if  we  use  Eq.  <G2),  wo  can  then  use  Eq,  (58)  to  express  P^q  as  a func- 
tion of  the  airplane  altitude  h rather  than  the  angle  0, 

For  feed  horn  at  (x^,  z^)  Instead  of  (xQ,  aQ),  we  get  for  the  power  pattern 
including  multipath 
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| VI'*  ■»  > 


JS# 


4*  PM1<6) 


FW^gOij)  -ikSj  r iCx 


-ML 


l5i  -‘S4 

dx  dy - — e * A,  e + A',"  v e 

a.  1 


rr  Fiwgin.)  -ika.  r i?.  -i?4l 

JJ  dxdy ±2 — e 1 Bj  e 1 + By  v e *| 

So 

rr  F((t1)g(n1)  r i€j  -i?4l 

JJ  dx  dy -*-g  ---  [C1  e 1 + cy  v e *J 


where 


8j  » [(x  - Xj)2  + y2  + (z  - Zj)2] l/2 


V (tti  + •»!  i)003  0 * 

Ay  - - ^ + yx  coa  - 02  + 6)  , 

Cl»-etao(«1  + 71!£)  , 

c'l"  * (“l  + fj)  ain  ^'2  “ °2  + 6>  » 


etc  . 


4.  SYUMKTiUC  Ufcl'LECTOH  AND  HUItN  MTKHNS 


(63) 


Let  ua  now  suppose  that  the  reflector  and  the  feed  hom  are  symmetric  in  the 
y- direction.  Tlrnt  ia,  for  Uie  reflector  surface  z ® f(x,  y)  we  have 


h 
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f(x,  -y)  = f(x,  y)  . 


Also,  for  the  horn  y-plane  pattern  we  assume 


g(-n0>  = g(nQ)  . 


Then,  we  write  £q,  (21)  as  (for  = 0 plane) 


(68) 


, ikR 
[2*R  e 0 

W 


-VA  f 

/ ~yJx) 


F(0jg(nJ  -ik(s  -xsin0-zco30) 
dy e 0 


{x  Aq(x,  y)  + y Bq(x,  y)  + z CQ(x,  y)} 


(69) 


?b  ya(x) 


/ 

xa  0 


F(0  )g(n  ) -ik(s  -xsinO-zcosO) 
dy— a— °-e  ° 


• {X  A0(x,  y)  + y BQ(x,  y)  + z CQ(x,  y) } , 

where  x^  and  are  the  coordinates  of  the  boundary  of  the  horn  along  the  x-axis 
and  yQ(x)  is  the  coordinate  of  the  boundary  along  the  y-axis,  as  shown  in  Figure  6. 


Figure  8.  Projection  of  z » f(x,  y)  onto  the  x-y  Plane 
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Now,  in  the  first  Integral  in  (69)  we  change  y to  -y  and  use  the  fact  that  for 
* = 0 (Note  thatao(-y)  = y).  PJ-y)  = -pQ(y)  and  yji-y)  = ^(y)] 


V-y ) = A0(y>  * 


Bo<-y> » -B<y>  . 


V-y^^y*  • 


30(-y)  = 3o(y)  . 


g[h0(-y)]  = g[n0(y>]  . 


This  gives 


2»R  e 


*0\  y»<x> 

-)h3  = 2 f d xj  dy  ^g.!!o!;lk(3o'X3y~2co3e) 

* ..  A O 


{*A  <x,y)  + 2c  (x.y)}  . 


Also,  from  Eq.  (62)  we  have  that 


A"'<-y)  ■ A^(y)  , 


B”*(-y)  » -J3”»(y)  , 


C”'(-y)  • C^'(y)  . 

Therefore,  from  Eq,  (26)  we  got  for  the  vacuum  pattern 
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where 


p = s - x sin  0 - z cos  0 . 

ro  o 


Also,  for  the  case  of  multipath  and  spherical  Earth  Eq.  (58)  becomes 
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The  resulta  Pj[0(h)]  and  PM1[0(h)] , which  we  would  get  for  the  case  when  the  feed 
horn  is  located  at  (Xj,  0,  Zj)  Instead  of  (xo,  0,  zq),  can  be  obtulned  from  Eqs.  (79) 
and  M)  by  replacing  iQ,  nQ,  VQ,  aQ,  AQ,  CQ,  A”',  C”' by  t)v  Pv  Sj,  A^ 

C A"'  ti'" 

Ul*  *1  * W * 

As  a check  on  the  validity  of  Eq,  (80),  we  can  consider  the  limit  when 
6 « tl'r,  = 0 (see  Figure  7),  In  this  case,  we  know  that  because  the  electric  field  is 
horizontally  polarized,  Uw  power  pattern  at  0>>  * -0  should  bo  zero.  In  tills  case, 
we  find  (for  aQ  — co) 

C*2  « 0"  * -0  , 

*3  * Ro  • 

*4  2 *1  • 

1 (81) 

AHI  * -A 
o o • 

C"»  * -C 
o o * 

v<  1 . 

So  tttat  upon  using  Uie  results  o ( Eq,  (81)  lit  (80)  gives 

pm0<8  . -d"  . -e2) . o , 

as  expected. 
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Figure  7.  Limiting  Case  of  Grazing  Incidence. 


5.  FURTHER  SIMPLIFICATION 

If  we  now  substitute  for  Aq,  Cq.  A^>,  C’"  etc.  in  Eqs.  (79)  and  (80),  we 
obtain  the  following  simplified  results: 


*V  (0) 


tf— ■ iv«!# 


*2P»n<®) 


$r—  * ly#|2  + v2iy-o")l2 


-2  v cos  -9g  + a - 0)  Be{y0)l*(-0")  elk^} 


where 


5>  ».w 


y«).^  <1*^  dy-i£~i!  (»0*vis) 


exp  {-1W«0  - a s in  9 - z cos  0) } , 


AR  = Rg  - Rq  = R1  + R2  - Rq  » 


Re  denotes  "real  part  of." 


6.  ERROR  IN  ALTITUDE  DUE  TO  MULTIPATH 


In  the  absence  of  multipath,  it  is  possible  to  obtain  information  on  an  airplane 
altitude  by  placing  two  horns  on  the  reflector  system  as  shown  in  Figure  1,  and 
then  operating  the  system  in  a monopulse  mode.  In  the  absence  of  multipath,  the 
voltage  V(h)  (see  Figure  1)  is  given  by 


V(h)  = iqiog  P^h)  - log  PQ(h>] 

rpi<h)1 

=ki°sLv5J  • 


where  K is  a constant,  Pj(h)  is  the  power  pattern  due  to  the  horn  at  (x^,  o,  z^), 
and  P0(h)  is  the  pattern  (see  Eq.  82)  due  to  the  horn  at  (xq,  0,  zq).  A typical 
curve  of  V versus  h (for  a given  airplane  range  r)  is  shown  in  Figure  On.  Now, 
because  the  measurement  of  V has  some  error  AV  in  it,  there  is  a corresponding 
error  in  our  calculation  of  tha  airplane  altitude  h.  This  is  indicated  pictorially  in 
Figure  8b. 

The  inclusion  of  multipath  compounds  our  problem  because  now 


vM(W 


K log 


LPM0<h>J  ' 


This  leads  to  an  additional  error  AV,,  given  by 


Figure  8a.  Plot  of  Output  Voltage  V Versus  Target  Altitude  h in  the  Absence 
of  Multipath,  Showing  How  Measurement  Error  AV  Leads  to  Altitude  Error 
Ah 


r a CONST. 


A typical  measurement  of  VM(h)  is  shown  in  Figure  9a  for  a fixed  target  range  r. 
The  envelope  of  AV^  is  also  indicated  on  Figure  9a.  Let  us  suppose  that  the  error 
for  a given  h,  say  hj,  is  AV^h^).  We  then  plot  this  error  on  Figure  8a  to  obtain 
the  result  shown  in  Figure  9b.  By  using  Figure  9b.  it  is  possible  to  calculate 
curves  of  the  error  Ah  (due  to  multipath)  in  altitude  versus  altitude  h for  a number 
of  different  target  ranges . A typical  result  is  shown  in  Figure  10a.  It  if  often 
convenient  to  plot  the  results  in  Figure  10a  versus  range  with  altitude  as  a param- 
eter; a typical  result  is  shown  in  Figure  10b. 

In  part  2 of  this  report,  we  will  employ  Eqs.  (82)  and  (83),  etc. , to  present 
curves  of  the  error  in  altitude  Ah  versus  h and  r for  specific  radar  configurations. 


ENVELOPE  OF  AVM(h) 


r = CONST 


